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o 

Q 

' We examine periodic, spherically symmetric, classical solutions of SU(2)- 

Higgs theory in four-dimensional Euclidean space. Classical perturbation 
theory is used to construct periodic time-dependent solutions in the neigh- 



borhood of the static sphaleron. The behavior of the action, as a function 
\ of period, changes character depending on the value of the Higgs mass. The 

' required pattern of bifurcations of solutions as a function of Higgs mass is 

examined, and implications for the temperature dependence of the baryon 
. number violation rate in the Standard Model are discussed. 

I. INTRODUCTION 

Ph. 

CD . 

Non-perturbative processes in quantum field theory such as tunneling, decay of 
metastable states, and anomalous particle production, are intimately connected with the 
existence and properties of non-trivial solutions of the associated classical field equations. 
In electroweak theory, baryon number violating processes are a consequence of topologi- 
cal transitions in which there is an order one change in the Chern-Simons number of the 
SU (2) gauge field. At zero temperature, such transitions are quantum tunneling events, and 
the rate of these transitions is directly related to the classical action of instanton solutions 
of SU(2) gauge theory Q. At sufficiently high temperatures,[] the dominant mechanism for 
baryon number violation involves classical thermally activated transitions over the potential 
energy barrier separating inequivalent vacuum states. The configuration characterizing the 
top of the barrier is the static sphaleron solution of S'f/(2)-Higgs theory ||]; the energy of 
this solution controls the thermally-activated transition rate 

As briefly reviewed below, when one lowers the temperature from the sphaleron dom- 
inated regime, the topological transition rate is related to the action of periodic classical 
solutions of the Euclidean field equations with a period (3 equal to the inverse temperature. 



But below the critical temperature (or cross-over) where "broken" electroweak symmetry is restored. 
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Very little detailed information is known, however, about the properties of periodic classical 
solutions in electroweak theory. 

In this paper, we examine periodic solutions in S77(2)-Higgs theory (which represents the 
bosonic sector of electroweak theory in the limit of vanishing weak mixing angle). We focus 
exclusively on classical solutions which are real in Euclidean space. Classical perturbation 
theory is used to construct periodic time-dependent solutions in the neighborhood of the 
static sphaleron. The variation in the period as one changes the amplitude of oscillation 
away from the sphaleron (or equivalently the turning point energy of the solution) is found 
to change sign depending on the value of the Higgs mass. Unlike the situation in many 
simpler models, there is no periodic classical solution which smoothly interpolates between 
instantons at low temperature (long period) and the sphaleron at high temperature (short 
period). We argue that, depending on the value of the Higgs mass, one or two different 
branches of periodic solutions must exist, connected at a bifurcation point. The topological 
transition rate (in the semi-classical approximation) must show a "kink" as the temperature 
varies between zero and the electroweak transition temperature. 

The non-existence of arbitrarily long period instanton-antiinstanton solutions in 577(2)- 
Higgs theory, and the consequent necessity for abrupt changes in the topological transition 
rate, has been previously noted [§,§J. This phenomenon also occurs in two-dimensional 0(3) 
non-linear sigma models with soft symmetry breaking terms, which mimic many features of 
577(2)-Higgs theory. Periodic classical solutions in softly-broken 0(3) sigma models have 
been studied in detail by several authors |5||||. We will argue that the behavior of periodic 
solutions found in the sigma model mimics the situation in S77(2)-Higgs theory when the 
Higgs mass is sufficiently small, but that for larger Higgs mass an additional bifurcation 
emerges. 

The plan of this paper is as follows. Section II briefly discusses classical solutions in 
a simple example of a one-dimensional periodic potential, and shows how bifurcations in 
solutions appear as one deforms the potential. This provides a useful analogy for the later 
discussion of changes in periodic i>77(2)-Higgs solutions as the Higgs mass is varied. Section 
III defines our notation for classical S77(2)-Higgs theory and its reduction to a 1+1 dimen- 
sional theory for spherically symmetric configurations, and summarizes the action of the 
relevant discrete symmetries. Periodic classical solutions resembling small, widely separated 
chains of instantons and antiinstantons are the subject of Section IV. Such solutions exist 
for periods small compared to the inverse mass scales of S77(2)-Higgs theory, and may be 
perturbatively constructed starting from superpositions of very small pure gauge instantons 
and antiinstantons. Section V contains a general discussion of classical perturbation theory 
for periodic solutions which are small deviations away from a static solution, and section 
VI discusses the practical issues which arise when applying this formalism to SU (2)-Higgs 
theory. Numerical results are the subject of Section VII, and the stability of periodic so- 
lutions close to the sphaleron is examined in Section VIII. The final section discusses the 
implications of these results, combined with the instanton-antiinstanton analysis, on the 
global structure of periodic classical solutions and on the temperature dependence of the 
(semi-classical) topological transition rate. 

Before continuing, we briefly review the relation between finite temperature transition 
rates and periodic Euclidean classical solutions. The partition function Z((3) for a quantum 
system at non-zero temperature (3~ l may be written as a Euclidean functional integral over 
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configurations which are periodic with period f3, 



Z{(3) = / [Vq\ exp(-S E [q]), (1.1) 

Jq(O)=q(0) 

where Se is the Euclidean action. For a semi-classical treatment of a weakly coupled theory, 
the classical action Se ~ ^/g 2 , where g is the small coupling constant, and the functional 
integral can be estimated by the method of steepest descent applied to the minima of Se- 
Non-perturbative transition rates may be shown to be related to saddle-point expansions 
about extrema of S E [q] possessing one negative mode |7|-|T0l. In particular, Euclidean peri- 



odic solutions with one negative mode determine the rate of quantum tunneling events at 
finite temperature. Generalizing the usage of Coleman J7|, we will generically refer to such 
classical solutions as "bounces" . To exponential accuracy, the number of tunneling events 
per unit time per unit volume, at a temperature scales as 

T ~ exp(-S E (l3)), (1.2) 

where Se{(3) is the Euclidean action of the bounce solution with period (3. 

II. A ONE-DIMENSIONAL EXAMPLE 

The behavior of solutions in a toy model of a one-dimensional periodic potential will 
provide a useful analogy.^ Consider the dynamics for a single degree of freedom defined by 
the Euclidean action 

S E [q] = Jdt Qg 2 + ^/(g)), (2.1) 

where the potential V(q) will be specified momentarily. Extrema of this action have a 
conserved energy E = —\q 2 + V(q) ) and it is a straightforward matter of integration to solve 
for periodic classical solutions. In particular, the period (3 is given by 

(3{E) = <f dq (2.2) 



where the integration takes place over a periodic trajectory with turning points qt p defined 
by E = V(qt v ). The Euclidean action of the classical solution with period (3 is 

S E ((3) = (3E(P) + jdq y/2[V(q)-E(j3)], ( 2 -3) 



where E((3) is obtained by inverting the equation for the period ( |2.2|) . 

If one varies the classical action S of periodic solutions, as a function of the period j3, 
the derivative dS/d/3 is minus the corresponding Hamiltonian (see, for example, Jill), f rom 
whence 



2 Similar discussion of much of this material may be found in Refs. [§,1). 
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(2.4) 



Consider the potential of a modified pendulum 



V(q) = ~f {[1- cos(g)] + A[l- cos(2g)]} , 



(2.5) 



with < A < 1/4. For A = 0, this reduces to a sinusoidal potential, while larger values 
of A flatten the maxima and sharpen the minima of the potential. Plots of the potential 
for A = 0, 1/8, and 1/4 are shown in figure [|. 
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FIG. 1. The modified pendulum potential V(q) (2.5), plotted for three values of the parameter A. 
For A = 0, it is the sinusoid shown as the solid curve. For A = 1/8, the maximum at q = 7r is flattened, and 
the minimum at q = sharpened, as shown by the dashed curve. This trend reaches its zenith at A = 1/4, 
where, as shown by the dotted curve, the quadratic term at the maximum vanishes. Larger values of A 
would cause the maximum to split so that V(q) would have two maxima per cycle. 

To understand the branching behavior of bounce solutions in this theory, one should 
examine the behavior of the potential V(q) near the sphaleron at q = it. The Taylor 
expansion of the potential about n begins 



V(n+x) = ^ 



;i-4A)^ + (l-16A)^-0(* 6 ) 



(2.6) 



In the inverted potential — V(q) in which Euclidean solutions evolve, one finds harmonic 
motion for small perturbations about the single sphaleron at q = tt for A < 1/4. For A > 1/4, 
as noted above, the sphaleron splits in two and the solution at q = tc becomes a local 
minimum of the action. 

For A < 1/16, the quartic term in the expansion ( |2.6| ) causes the inverted potential to 
soften from what the simple harmonic term would give. Oscillations thus increase in period 
as their amplitude is increased, or alternatively, as the turning point energy is lowered from 
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the sphaleron energy. This trend continues all the way to vanishing turning point energies, 
where, at long periods, the bounce approaches a periodic kink plus antikink. 

For A > 1/16, however, one sees from the expansion ( |2.6| ) that the quartic term agrees 
in sign with the quadratic term in the potential, and thus causes the potential to be steeper 
than what the simple harmonic term would give. This implies that the bounce, at least 
initially, decreases in period as its amplitude is increased, or as the turning point energy is 
lowered from the sphaleron energy. At some critical turning point energy, this trend must 
reverse in order to connect the bounce solution to the kink-antikink at very long periods. 

As noted above, this model is straightforward to integrate numerically. The results are 
shown in figure |2|, for three different values of the parameter A. The plot shows the action 
of Euclidean classical solutions of the theory as functions of the period. As indicated by the 
derivative of the action (|2.4|), the slope of a curve on such a plot is the turning point energy 
of the corresponding classical solution. The dotted line is the sphaleron solution, which, as 
a static configuration, has the same energy at any value of the period. 

At the critical value of A = 1/16, the bounce still increases monotonically in period 
as its turning point energy is decreased, giving the single curve of solutions starting from 
the sphaleron at point Q. Above this value of A, the bounce initially decreases in period 
as its turning point energy is decreased, then reaches a bifurcation point, after which the 
bounce period begins to increase again as the turning point energy is further decreased. 
For A = 3/32, the bounce begins at the sphaleron at point R, then proceeds to the bifurcation 
point X as the turning point energy is lowered, after which further decreases in turning point 
energy take one towards longer periods. As A is increased, the length of the branch RX 
increases, becoming the branch SY at A = 1/8. 
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F IG. 2. Curves of action vs. period for periodic "bounce" solutions in the modified pendulum poten- 
tial (2.5), in units where E$ — 1. The slope of a curve is the turning point energy. The straight dotted 
line represents the static sphaleron. The curve merging with the sphaleron at point Q represents bounces 
with the critical value of the parameter A = 1/16. Above this value of A, the curves split at bifurcation 
points. At A = 3/32, for example, the bounces pass from point 7? at the sphaleron, towards shorter periods, 
until reaching point X, beyond which further decreases in the turning point energy take one towards longer 
periods. The branch RX grows longer for larger values of A, becoming the branch SY for A = 1/8. 
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The (linearized) stability of a classical solution is determined by the number of eigen- 
modes of the curvature of the action S 2 S (evaluated at the solution of interest) with negative 
eigenvalues. Bifurcations, such as X or Y in figure |2|, correspond to points in the space of 
solutions where a curvature eigenvalue passes through zero. When different branches of so- 
lutions merge at a bifurcation such as Y, the branch of solutions with larger action has one 
more negative mode than the lower branch. In the case of the pendulum, the negative mode 
of the upper branch can be identified with a deformation which moves the turning points 
further in or out, while holding the period fixed. Such a deformation causes no first order 
change in action, since the momentum vanishes at the turning point. But to second order, 
such a deformation can cause the action to change. The eigenvalue of this mode passes 
through zero at the bifurcation point, and becomes positive on the lower branch of solutions 
emerging from the bifurcation. 

Points Q, R, or S represent bifurcations where a branch of periodic bounce solutions 
emerges from the static sphaleron. Each point on the branch SY (or RX) represents a one- 
parameter family of solutions related by time translations, i.e., x(t— to). The time deriva- 
tive x(t— to) is a zero mode in the small fluctuation spectrum and represents a deformation 
which is an infinitesimal shift in to- The sphaleron, being static, has no such zero mode. 
At periods longer than that of the bifurcation S (or R) the sphaleron has three negative 
modes: one is the static unstable mode, corresponding to a deformation x — > x + e with time- 
independent e, while the other two are periodic oscillations, 5x oc sin(27rt//5) or cos(27rt//3). 
At the bifurcation point, the two oscillating negative modes become zero modes, and at 
shorter periods, positive modes of the sphaleron. On the branch SY (or RX) of bounces 
emerging from the sphaleron, one of the zero modes remains the time-translation zero mode 
of the bounce, while the other becomes a negative mode — the same negative mode which 
later passes through zero at the bifurcation Y (or X) at the other end of this branch of so- 
lutions. The bounce solutions also have one additional negative mode resembling the static 
unstable mode of the sphaleron. 

This connection between bifurcations and stability of solutions will be important for the 
discussion in section IX. Several figures illustrating the "topography" of the action near 
bifurcations of the ,S£7(2)-Higgs sphaleron, and a more detailed discussion of the negative 
and zero modes of the various solutions, appear in Section VIII. 

At any temperature, the solution which characterizes the most probable barrier crossings 
(and whose action controls the rate of barrier crossing) is the bounce, or sphaleron, with the 
smallest action for period f3 = 1/T. As is clear from figure |2], this means that the branches, 
such as RX and SY, which stay above the action of the sphaleron, do not control the rate 
of barrier crossings at any temperature. Instead, for A > 1/16, there is an abrupt change 
from sphaleron-dominated barrier crossing to bounce-dominated crossing at the tempera- 
ture where the curve of bounces crosses the sphaleron line. Physically, this is a transition 
between classical thermally-activated transitions over the barrier, for temperatures above 
the crossover, and quantum tunneling with a most-probable energy E < E$ for tempera- 
tures below the crossover In contrast, for A < 1/16, as one increases the temperature from 



3 In the leading WKB approximation, the transition rate is discontinuous at the crossover. This discontinuity 
is smoothed into a narrow transition region of width 0(g 2 ) in the exact transition rate. 
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zero the most probable energy for tunneling grows until it reaches E$, at which point the 
turning points of the most probable tunneling path merge at the top of the barrier (and 
the WKB transition rate smoothly interpolates between quantum tunneling and classical 
thermal activation ||). 



III. THE S77(2)-HIGGS MODEL WITH SPHERICAL SYMMETRY 

We now turn to the SU(2) gauge theory in 3+1 dimensions with a single Higgs scalar in 
the fundamental representation. This model represents the bosonic sector of the Standard 
Model of the weak interactions in the limit of small weak mixing angle. The action may be 
written in the formf] 



S = — I d l .r 



1 M 2 

-Tr(F^F^) + {D^YD^ + tttJ-^ - 2M 2 W ) 

i qM w 



(3.1) 



The scalar field <£> is an SU(2) doublet with covariant derivative = (d^ + A^), where = 
A a ^r a /2i. The gauge field strength is the commutator F M „ = [D^,D U }. 

The action ( |3.1|) has an explicit SU{2) L gauge symmetry, represented by SU{2) matrices 
acting on the Higgs doublet $ (or its conjugate $ c = — ir 2 $*) from the left. It also has a 
custodial SU(2)ji global symmetry, given by SU(2) matrices multiplying the matrix ($, $ c ) 
from the right. 

Spherically symmetric field configurations are those for which the effect of a rotation 
can be undone by a gauge transformation combined with a custodial SU{2)r transforma- 
tion. Imposing spherical symmetry reduces the four- dimensional theory to a two-dimensional 
theory [|T2| , p~3|] , which can be parameterized in terms of six real fields a, /3, /i, z/, a , and 
a\. In terms of these two-dimensional fields, one may write the original fields of the four- 
dimensional theory as 

A (r, t) = a f • r /2i , (3.2a) 
1 



-(a— 1) t x f H — (3 (r — (f • r)r) + ai(r • r) f 



(3.2b) 

$(r,t) = (fx + iur -f)£. (3.2c) 



Here f = f/r is a unit radial vector, and £ is a unit doublet which can be rotated arbitrarily 
by combined global SU{2) L and custodial SU{2) R transformations. It is helpful to define 
complex linear combinations x = a + i/3 arid <fi = \i + iv '. The field \ then represents the 
spherically symmetric degrees of freedom of the tangential gauge fields, and those of the 
Higgs field. 

After imposing spherical symmetry in this form, the two-dimensional theory which re- 
sults has a U(l) gauge symmetry, which is the subgroup which remains of the original 



4 For convenience, we have rescaled the gauge and Higgs fields so that the action has an overall 1/g 2 - 
The conventional Higgs vacuum expectation value is v = 2Mw/g, and the usual quartic coupling is A = 
gg 2 Mfj J 'Myy . As usual, the weak fine-structure constant is aw — 3 2 /47r. In Minkowski space (with a 
space- like metric), the usual action is minus (B.l). 
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SU(2)l group. As elements of the original SU(2) group, these U(l) gauge transformations 
are given by Q(f,t) = exp(iu(r,t)T ■ f/2). Under these Abelian gauge transformations, 



(f> — > exp(iuj/2)(f), 

It is thus natural to define two-dimensional covariant derivatives as 



(3.3a) 
(3.3b) 
(3.3c) 



D^x = (<9 M - idp) x, 
= (dp - ictp/2) . 



(3.4a) 
(3.4b) 



Here, the reduced space-time indices \i and v take on the values and 1, corresponding to 
the time t and distance from the origin r, respectively. We also define a (1 + 1) -dimensional 
U(l) gauge field strength as f^ v = d^a u — d v a^. 

Substituting these definitions into the action ( |3.1| ) and integrating out the angular coor- 
dinates gives for the two-dimensional action 



S 



4:71 
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dr dt 



4 



r 2 f, u r+\D X \ 2 + r 2 \D^ 



| X | 2 +l)-Re(x 



M H 2 



2 -2M 2 v f 



(3.5) 



1 as r 



0, and 



Note that finite action configurations must satisfy 
y/2Mw, \x\ ~> 1 with x*0 2 re& l an d positive as r — > oo. 

The transformations of the four- and two-dimensional fields under the discrete symme- 
tries of the S7/(2)-Higgs model are summarized in Table |I[ When acting on spherically 
symmetric configurations, four-dimensional charge conjugation can be undone by a global 
SU(2)l rotation combined with a custodial SU(2)r rotation, and therefore it has no effect 
on the fields of the two-dimensional reduced theory. Instead, U(l) charge conjugation in the 
two-dimensional theory is produced by a parity transformation in four dimensions. 

The sphaleron is a static field configuration for which the action is stationary. It may 
be found, in one choice of gauge, by imposing £7(1) charge conjugation invariance on the 
spherically symmetric fields defined above. Referring to the action of parity in table |, one 
sees that U(l) charge conjugation invariance implies that a M = 0, and </> and x are purely 
real. There are thus only two real fields, a and \l. The sphaleron is the lowest energy 
configuration for which a changes from —1 to 1 as r goes from zero to infinity. The static 
energy to be minimized can be written, by simplifying the above action, as 



E 



47T 
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dr 



(d r a) 2 +r 2 (d rf x) 2 + ^(« 2 -l) 2 + ^(a-l)V + 



M 2 H 



-2M 2 W 



(3.6) 



Making an appropriate substitution on the radial variable r to map the positive real line onto 
a finite interval, and then discretizing the above expression, one can readily solve numerically 
for the sphaleron fields a sp h and /i sp h [@,[L3|,|l4[ • 



S 



Field 


C 


P 


T 


$(f,t) 
A{r,t) 
Ao(f,t) 


$*(f,t) 
A*{r,t) 
A*(r,t) 


$(-r,t) 
-A(-f,t) 
A (-f,t) 


$(f, -i) 
A(f, -t) 
-A (r,-t) 


<f>{r,t) 
X{r,t) 
ai{r,t) 
a (r,t) 


4(r,t) 
X(r,t) 
ai{r,t) 
a (r,t) 


<f>*(r,t) 
X*(r,t) 
-ai(r,t) 
-a (r,t) 


</>(r, -t) 

-*) 

ai(r, -t) 
-a (r, -£) 



TABLE I. The transformations of the various four- and two-dimensional fields under the action of the 
discrete symmetries of the SU (2)-Higgs model. C is charge conjugation, P is parity, and T is time rever- 
sal. Charge conjugation C can be compensated by a global SU(2)l rotation combined with a custodial 
SU(2)r rotation, and so becomes an identity transformation on the two-dimensional fields. Parity P pro- 
duces the effect of U(l) charge conjugation on the two-dimensional fields. 



IV. INSTANTON-ANTIINSTANTON SOLUTIONS 



The 5'?7(2)-Higgs model discussed in the preceding section does not have exact instanton 
solutions, that is, finite action solutions on unbounded Euclidean space. Derrick's theorem 
(see for example fl5]) shows that a reduction in the four-dimensional length scale of any 
purported instanton will decrease its action. There are, however, approximate solutions 
resembling instantons of the pure SU(2) theory with characteristic sizes much smaller than 
the inverse mass scales of the theory, namely M^ 1 or Mjj 1 [|TJ. The action of these con- 
figurations decreases as the instanton shrinks (asymptotically approaching the pure gauge 
value of 87r 2 /g 2 ), but in all other directions the action is minimized. At zero temperature, 
these small instantons dominate topological transitions in SU (2)-Higgs theory; quantum 
fluctuations stabilize the size of the relevant instantons 0. 

Exact classical solutions resembling a periodic chain of small instantons and antiinstan- 
tons also exist in SU (2)-Higgs theory^ Such a configuration is sketched in Figure |3]. In 
order for this configuration to be a classical solution, the attraction between the instan- 
tons and antiinstantons must exactly balance the tendency of the individual instantons or 
antiinstantons to collapse. The resulting solutions may be constructed iteratively starting 
with a superposition of SU (2) pure gauge instantons, provided the instanton sizes are small 
compared to their separation, and the separation is small compared to the inverse mass 
scales Mh 1 or M T 7 X 



5 Much of the following material may be found in Ref. |Tc|| . 
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FIG. 3. A sketch of periodic instanton-antiinstanton configurations. The r axis represents the spatial 
distance from the origin, while the r axis is imaginary time. The circles / are instantons, of size p, while the 
circles / are antiinstantons, of size p. The distance from an instanton to the nearest antiinstanton is t, while 
the period is (3. To construct such a solution iteratively, p and p must be small compared to t, (3, or /3—t, 
which in turn must be much shorter than the inverse mass scales of the theory, labeled here as M . 

To leading order, the fields of an instanton centered at the origin in singular gauge are JIJ 

M*) = x2(p 2 2 ^ x2) %u^(j a M [1 + 0(M w x)), 

= V2Mwx^ + 0{Mwx) + 0{Mhx)1 ( 41 ) 
\Jx l + p 2 

where £ is an arbitrary complex unit doublet, and 77* denotes 't Hooft's 77 symbol The 
action of this field configuration, to leading order in the small parameter Mp (where M 
denotes either My/ or Mh) is 

Si = % + ^-M^p 2 + 0(M*p*) . (4.2) 
9 9 

The 0{M^ v p 2 ) correction to the zero-order instanton action is due to the derivative term for 
the Higgs field, and is readily calculated as a surface integral QTJ] . The neglected C(M 4 p 4 ) 
contribution comes from the Higgs potential. 

Consider placing such an instanton and antiinstanton along the time axis a distance 
t apart. Since the fields approach their vacuum values for large distances t ^> p, one 
can construct an approximate solution to the equations of motion by linear superposition. 
Furthermore, if p,p t M™ 1 , the leading interaction term between an instanton and an 
antiinstanton has the dipole-dipole form of the pure gauge theory flEfl . It is 



ASn 
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(4.3) 



when the group orientations of the instanton and antiinstanton are aligned to give the 
maximally attractive interaction. This dipole interaction arises from the pure gauge field 
part of the action, and may be calculated as a simple surface integral. The interaction term 
for two instantons, or two antiinstantons, vanishes to this order. 
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For a periodic chain of alternating instantons and antiinstantons, such as the one illus- 
trated in figure there is a dipole interaction between each instanton and antiinstanton. 
Summing over all of these interactions yields, for the action per period of the periodic 
instanton-antiinstanton, 



Sn(p, P, t; p) 



g 2 

4tt2 
9 2 



4 + 



A + M 2 w {p 2 +p 2 )- 



H n=—<y 
^ P 2 p 2 ( 



24 



3-2sin 2 (7^^//5) , 
sm\nt/[3) 



(4.4) 



neglecting corrections of order p a //3 6 , M 4 p 4 , and Mwp 4 /j3 3 . This expression is stationary 
with respect to variations of t when t = (3/2. Requiring the derivatives with respect to p 
and p to vanish then fixes 



P = P = 

so that the classical action per period is 



V2M W (3 2 

471-2 



r 



M 4 tf 4 

16vr 2 + ^f- + 0(M^/? 5 ) 



(4.5) 



(4.6) 



This periodic solution has two negative modes. One negative mode, due to the attraction 
between the instanton and antiinstanton, is obvious because the action (|4.4| ) is maximized 
with respect to t. The other negative mode is a symmetric perturbation 5p = Sp in instanton 
and antiinstanton sizes. This can be easily verified by noting that the mixed second derivative 



d 2 S 



dp dp 



16tt 2 M 2 



w 



extremum 



(4.7) 



is the only second derivative of S involving p or p which doesn't vanish at the extremum of S. 
The right hand side of the mixed second derivative (|4.7| ) is thus a negative eigenvalue of the 
curvature of S, with an eigenvector satisfying 5p = 5p (and St = 0). This second negative 
mode is a signature of the delicate balancing required to construct a genuine periodic solution 
out of widely separated instantons, which themselves are not solutions. 

The periodic instanton-antiinstanton solutions are time-reversal even, and thus have 
turning points, where all of the time-reversal odd field components (and velocities) vanish. 
The turning points are located halfway between the instanton and antiinstanton (at the ori- 
gin in figure ^|), as is clear by observing that the time-reversed fields of an instanton are those 
of an antiinstanton, and vice versa. The turning point energy of these solutions vanishes 
as the period goes to zero. Both negative modes of the periodic instanton-antiinstanton 
solutions are parity odd and time-reversal even. 

V. THERMAL "BOUNCES" NEAR THE SPHALERON 



The 5'?7(2)-Higgs action (|3.1|) is time-reversal invariant. Solutions to the equations of 
motion must either be time- reversal invariant, or come in time-reversed pairs. The periodic 
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bounces we are interested in are time-reversal invariant solutions. This means that, at some 
time during the period, all time-reversal odd quantities must vanish. These times correspond 
to classical turning points, at which the kinetic terms in the Lagrangian vanish, leaving just 
the potential energy. 

At turning-point energies just below the sphaleron energy, the bounces will be small 
fluctuations away from the sphaleron. Therefore, they may be studied using classical time- 
dependent perturbation theory. Given a time-independent action functional S[<p] of fields 
generically denoted as 0, and a static solution O (the sphaleron) which satisfies the equation 
of motion 5S , [0 O ] = 0, consider adding a small periodic fluctuation 50, of period (3, to O . 
The classical equation of motion can be expanded in powers of the small fluctuation 5(f): 

= 5S [0] = 5S [0 O + 50] 

= S (2) " ¥ + ±S® ■ (50) 2 + i S< 4 > • (50) 3 , (5.1) 

where the inner product implies summation over internal indices as well as integration over 
space-time, and we have defined a shorthand notation for variational derivatives of the action 

= 8 n S[(f) \. (5.2) 

The leading order piece of the equation of motion flS.lp is 

O = S {2) -50 + e>((50) 2 ). (5.3) 

This implies that the curvature has (at least) one eigenvalue of magnitude (9(50); that 
is, which becomes zero in the limit of vanishing fluctuation amplitude. This zero crossing of 
a curvature eigenvalue occurs at, and defines, the critical period (3q. The presence of a zero 
mode at the bifurcation point was pointed out in the context of the pendulum model; here 
we see its necessity from perturbation theory. We will label this particular eigenvalue A(/3), 
explicitly indicating its dependence on the period of the fluctuation, so that 

A(/3) = A'(/3 ) 09 - AO + 0(((3-(3 ) 2 ) . (5.4) 

Let rj(P) denote the corresponding eigenmode, so that 

. = A(/3) ntf) . (5.5) 

To leading order, 50 consists of an undetermined (small) amount of the zero mode 77(A))- 

To solve for 50 iteratively, beyond this linearized order, we rewrite the equation of 
motion using the projection operator onto the eigenmode r], denoted P v , to separate 
the small eigenvalue A when writing the inverse curvature = PJX + S® (1-P„). 

This gives 

50 = -^)- 1 -(i5( 3 )-(50) 2 + i^)-(50) 3 ) 

= ev~ S {2r \l - P v ) ■ Q ■ (50) 2 + I fi<<) ■ (50) 3 ) , (5.6) 
where the second line of this expansion ( |5.6|) follows from the definition 
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jV ■ Qs< 3 > ■ (5^ + i ■ (50) 3 ) . 



(5.7) 



(For convenience, we have assumed that the eigenmode i](/3) is normalized to one.) De- 
fined in this way, e is of magnitude 0(5<f)), and will be used as a small parameter in the 
expansion (|5.6|) . 

The last line of the expansion (5.6) defines a recursive recipe for computing 5(f) as a 
power series in the small parameter e. Each place 5(f) occurs in the last line, one may 
substitute an entire copy of the last line, stopping at whatever power of e is appropriate. 
This expansion is most conveniently written as a series of tree graphs, in which n-point 
vertices (for n = 3 and 4) correspond to — S^ n \ leaves ( — •) emerging from any vertex 
represent factors of erj, and internal lines correspond to the propag ator \l-P n )- Each 
graph is to be divided by the usual symmetry factor, given by the number of permutations 
of leaves or vertices which leave the graph unchanged. 



Carrying out the iterative expansion Q5.6| ) explicitly to 0(e 3 ) gives us 



+ 



I. 



+ 



0(6*) 



(5.8) 



The graphical expansion Q5.8D shows the beginning of the perturbative series for the field 
fluctuation 5<p in the small parameter e which may be used to characterize its size. But 
the parameter e is not independent of the oscillation period (3. The defining equation ( |5.7|) 
relates the value of e to the period j3, and may be used to generate a perturbative expansion 
of the period (3 in powers of e. To see how this goes, expand equation f[7\ using the same 
graphical expansion for 5(f), 



e 2 A 



K + 4 >"< 



+ 



0(e 



(5.9) 



In fact, all diagrams in the expansion ( |5.9|) with odd powers of e (and hence rj) vanish. This 
is because these diagrams imply integration over all time variables, and rj has a definite 
frequency as will be seen explicitly in the following section. To fully expand the left- 
hand side of the A equation ( |5.9|) , use the expansion ( |5.4|) for the eigenvalue A(/3) and then 
expand f3— (3q in powers of e. Since, as noted above, only even powers of e occur on the 
right-hand side of Eq. (|5.9p , the same must be true in the expansion of the period (3, 

(5.10) 



2!'" " 1 3! 

Substitute the period expansion ( |5.10| ) into the eigenvalue expansion ( |5.4| ), substitute the 
result into the A equation ( |5.9|) , and match corresponding powers of e. The leading order 
result is 



e 4 A'(/5 )/3«=4^--<] 



or more explicitly, 



A'(/3 ) 



1 



(5.11) 



(5.12) 



where all quantities may now be evaluated at the critical period (3q. 
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VI. COMPUTATION 



We will evaluate equation ( 5.12| ) for perturbations about the sphaleron in S*?7(2)-Higgs 



theory, in order to find the leading-order dependence of the period of sphaleron oscillations on 
the amplitude e of the fluctuating fields. We will also evaluate the leading order dependence 
of the turning point energy on the amplitude e in order to re-express dependence on e as 
dependence on the turning point energy. 

To carry this out, one must discretize the radial dependence in the action ( |3.5| ), solve for 
the static sphaleron, determine the critical period /3 and the leading order fluctuation rj(P ) 
by finding the zero-crossing eigenmode of the curvature evaluated at the sphaleron, and 
then assemble the various pieces needed to compute ( |5.12|) . 

Since we are looking for periodic solutions, the time dependence of each of the fields 
appearing in the field equation ( |5.8| ) may be expanded in a Fourier series. Because is a 
time-independent linear operator, only the vertices and generate couplings between 
different harmonics. For the calculation of f3^\ only field components at frequencies 0, fig 1 , 
and 2@q 1 are needed. 

Fixing a gauge is required, so as to remove zero modes of the curvature due to the residual 
U(l) gauge invariance. There are two obviously reasonable choices: radial gauge (ai = 0), 
or temporal gauge (a = 0). We chose to work in radial gauge, in order to avoid having to 
introduce link variables in the one-dimensional radial lattice. 

Because the action, and the sphaleron in our choice of gauge, are parity (or U(l) charge 
conjugation) even, it is useful to separate the fluctuating fields into parity even and odd 
components. The first order fluctuation rj is entirely odd. Each diagram in the field equa- 
tion (|5.8| ) then produces a contribution which is even if the number of leaves is even, and 
odd if the number of leaves is odd. The curvature is block-diagonal, and only the parity 



even block is needed for ( |5.12j ). Since the zero mode r]((3 ) is parity odd, the projection 
operator (1— P v ) has no effect on the parity even block of the propagator. 

Table [TI| shows the discrete symmetries and leading order perturbative expansions of the 
real fields appearing with this choice of gauge. Note that all time derivatives appearing in 
the action, or its variations, may be trivially computed analytically. 

The spatial derivatives in the action ( |3.5|) were discretized using a uniform sampling of 
the transformed radial variable |0| 



In 



1 + M H r 
1 + M w r 



\n(M H /M w ), (6.1) 



which maps the semi-infinite line < r < oo onto the unit interval < s < 1. Once the 
action is discretized, the sphaleron field, and the variations of the action at the sphaleron, 
can be computed as functions of the period. 



6 Temporal gauge has a small advantage of allowing one to reduce the Lagrangian to the separable 
form (d t <fi) 2 + V((f>), if 4> are suitably rescaled field variables. This simplifies the calculation of the critical 

period /3o, as — (27r//3o) 2 will simply equal the negative eigenvalue of the static curvature V^((f>). However, 
this choice has the cost of introducing the link variables into the discretization of all spatial derivatives 
on the radial lattice, and requires an otherwise unnecessary similarity transformation on the field variables 
to obtain the separable form. 
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Field 


P 


T 


Leading Order Terms 


ct(r, t) 


+ 


+ 


"spilth 


+ e 2 


«(°)(r) + a (2 \r) cos(2ut) 




fi{r,t) 


+ 


+ 


/isph(r 


) + e 2 


/i(°)(r) + /i^(r) cos(2cut) 




/3(r,t) 




+ 




e 


(3 {l \r) cos(ut) 


u(r,t) 




+ 




e i/^(r) cos(cut) 


a (r,t) 








ea^\r) sin(o;t) 



TABLE II. Symmetry properties, and leading order terms, of the field components for small oscillations 
about the sphaleron. P is parity (or equivalently U(l) charge conjugation), T is time reversal, and a + or — 
indicates whether the field is even or odd, respectively, under the symmetry. The static sphaleron fields are 
labeled a S ph and /i sp h- The fundamental angular frequency u) is 2tt/Pq, where (3$ is the critical period. The 
zero eigenmode ?y(/3o) is composed of the leading order terms for the parity odd fields /3, u, and a$. Note 
that only even powers of e occur in the expansions for parity even fields, and only odd powers of e appear 
in parity odd fields. 

In particular, the curvature of the action is a matrix which, if inverted at the critical 
period, with the zero mode projected out, would give the propagator. In practice, it is 
computationally inefficient to actually calculate the matrix inverse. Since the action is local, 
discretizing the radial derivatives as nearest-neighbor differences gives the curvature matrix a 
band diagonal structure. Band diagonal algorithms for finding eigenvalues and solving linear 
equations have a computational cost which scales only linearly with the matrix dimension. 



An algorithm for band diagonal matrices [17] was used to find the eigenvalues of the curvature 



matrix, including A(/3). Another band diagonal algorithm was used to solve the linear 



equations ■ z = £, with £ = ■ rf. The same linear band algorithm was employed to 
find the sphaleron fields /i sp h and by Newton's method iterations, and to generate the 
eigenvector rj(/3) by inverse iteration |TS| . 

As noted above, the variation (3^ of period with amplitude ( |5.12| ) may be converted into 
the variation in the period of the bounce with respect to the turning point energy, at the 
sphaleron. Explicitly, 



dE 



(6.2) 



Since the kinetic terms of the Lagrangian vanish at the turning points, the variation of the 
turning point energy with amplitude dE/de 2 is computed as the inner product 



dE 
Ik 2 



Wo) 



00 



ifcp(A0 • s® ■ ifo,(A0 



(6.3) 



where rj tp (p ) are static fields equal to the turning point values of rj(j3o). Similarly, the 
derivative with respect to period of the eigenvalue X'(Po) appearing in the (3^ equation 
( |5.12| ) is computed as the inner product 



A' (A,) = 77(A)) • • v((3o) ■ 



(6.4) 
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VII. NUMERICAL RESULTS 



Figure ^ shows the dependence of the critical period (3q of oscillations on the Higgs mass 
Mn. The critical period was found by adjusting an arbitrarily large period j3 down until the 
eigenvalue of the oscillating negative mode of the curvature matrix crossed zero. 

Figure |5| shows the action of the sphaleron for one critical period /?o as a function of the 
Higgs mass Mh- The sphaleron action was calculated by simply multiplying the sphaleron 
energy by the critical period (3q. This plot is important for understanding the connection 
between the sphaleron oscillations and the periodic instanton-antiinstanton solutions which 
occur at very short periods. The minimum action of these small instanton-antiinstanton 
solutions over one period, which is An /aw, is shown as the dashed horizontal line on the 
graph. The sphaleron action crosses this minimal instanton-antiinstanton action when Mh ~ 
6.665 M w . 

6.0 -| 



5.0 - 



(3 M W 4.0 - 



3.0 - 




0.0 2.0 4.0 6.0 8.0 10.0 

Mh /My/ 

FIG. 4. The critical period of sphaleron oscillations /3o vs. the Higgs mass Mh- 
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18.0 i 



16.0 



PoM^aw 14.0 



4k 

12.0 



10.0 H 1 1 1 — i 1 1 

0.0 2.0 4.0 6.0 6 665 8.0 10.0 

M H /M W 

FIG. 5. The action of the sphaleron over a single critical period vs. the Higgs mass Mjj. The dashed 
horizontal line shows the minimum action of an instanton-antiinstanton solution. 

Figure ^ shows the derivative of the period with respect to the turning point energy, 
d/3/dE, as a function of the Higgs mass. At low Higgs mass, the derivative is positive, 
meaning that the oscillations move towards shorter periods as the turning point energy is 
decreased. At a Higgs mass of 3.091 Mw, the derivative crosses zero and thereafter becomes 
negative, meaning that for Higgs masses larger than this value oscillations move toward 
periods longer than f3$ as the turning point energy is initially decreased from the sphaleron 
energy. 

1.0 -i 




-1.0 



-1.5 



0.0 



-I 1 1 1— 

2.0 3 - 091 4.0 6.0 

M H /M W 



8.0 



10.0 



FIG. 6. The derivative of the period of small sphaleron oscillations with respect to the turning point 
energy vs. the Higgs mass M H . 
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VIII. STABILITY 



Unstable modes of the bounce solution are closely related to those of the sphaleron. 
The sphaleron has a single static negative mode For periods shorter than the critical 
period /3 , this is its only negative mode. For periods (3 longer than /3 , the sphaleron has 
two additional negative modes with frequency uj = 2tt//3, given by the eigenmode r)(/3) and 
its time derivative d t T](d). 

Figure [7] illustrates the landscape of the action near the sphaleron and nearby bounce 
solutions, in the subspace spanned by 5(f) and its time derivative d t 8(f>. The action is repre- 
sented by the height of the surface in each graph. The periodic time translation symmetry is 
represented as the rotational symmetry around the vertical axis. Classical solutions appear 
as peaks, cups, or saddle points. The sphaleron sits on the symmetry axis, because it is 
time translation invariant. The direction corresponding to the static negative mode of the 
sphaleron is not shown. 

As shown in figure [Tj, when (3 < j3 , bounce solutions appear as a circular ridge around 
the sphaleron and have two negative modes (one of which is a small perturbation of the 
static negative mode of the sphaleron, and is suppressed in figure |7J), plus a time-translation 
zero mode. When (3 > (3q, bounce solutions appear as a circular valley around the sphaleron 
and have only one near-static negative mode plus the time-translation zero mode. Both 
the near-static negative mode of the bounce, and the additional negative mode for (3 < /3 
(which is a small perturbation of f]{f3o)) are parity odd and time-reversal even. 



(3<(3 



(3>(3 



dE 



(Po) > 





9(3 (3) 



< 





Ks = 1 



Kb = 1 



FIG. 7. Topography of the action near the sphaleron and nearby bounce solutions. The height of the 
surface in each graph represents the action. The extremum of the action on the axis of rotational symmetry 
is the sphaleron. The other extrema are the periodic bounce solutions, and have a time translation zero 
mode represented by rotations about the symmetry axis, ks denotes the number of unstable modes of the 
sphaleron, while Kg is the number of unstable modes of the bounce. The action has one additional unstable 
direction, not shown in the figure, corresponding to the static unstable mode of the sphaleron. 
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IX. DISCUSSION 



The results shown in figures |5] and [(], when taken together with the analytic results on 
short-period periodic instantons, allow one to suggest a picture of the branches of periodic, 
spherically symmetric, Euclidean classical solutions of this theory, which will be illustrated 
on the following plots of action versus period. On such plots, periodic classical solutions 
trace out curves, the slope of which at any point is equal to the conserved energy of the 
solution. The sphaleron, as a static solution, appears as a straight line through the origin, 
with a slope given by the energy of the sphaleron. Bounces with turning points have a slope 
equal to their turning-point energy, and thus must always have a positive slope. Solutions 
may merge (at bifurcation points where curvature eigenvalues pass through zero) when the 
curves representing the solutions meet at tangent points on the plot. The curves must be 
tangent for the solutions to have the same conserved energy. 



Action 




Period 



c) M H > 6.665 M w 

FIG. 8. Euclidean classical solutions of the ST/(2)-Higgs model, sketched on plots of action versus period, 
for three different values of the Higgs mass. On these plots, slope is equal to turning-point energy. The plots 
are drawn for purposes of illustration, and are not to scale. The line through the origin is the sphaleron. 
The curve of bounces C merges with the sphaleron at the point Q, which is at the critical period /3rj- For 
periods very short compared with the inverse mass scales (M^ 1 , M^ 1 ) of the theory, there exist periodic 
instanton-antiinstanton solutions / with near-zero turning point energy. 

Figure § summarizes the perturbative results for the bounces and periodic instanton- 
antiinstantons found in this paper, for three different values of the Higgs mass. These 
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schematic plots of action versus period are not drawn to scale. The curve of bounces C 
merges with the sphaleron line at point Q, which is at the critical period f3$. For periods very 
short compared with the inverse mass scales (Mjy , M w v ) of the theory, there exist periodic 
instanton-antiinstanton solutions / with near-zero turning point energy. For brevity, we will 
henceforth call these, in an abuse of language, periodic instantons. 

For light Higgs (less than 3.091 Mw), figure ||(a) shows that the bounces C decrease in 
period as their turning point energy is decreased from the sphaleron energy. This is clear 
from figure ^ for small Higgs mass. On the same plot, the action at point Q is larger than 
the minimum action of the periodic instantons /. This one can read from figure |5] in the 
limit of small Higgs mass. 

As noted in figure H, above a Higgs mass of 3.091 Mw, the bounces C increase in period as 
their turning point energy is lowered from the sphaleron energy. This situation is illustrated 
in figure ||(b). 

From figure [| we note that as the Higgs mass is increased, the action at point Q decreases, 
because the period (3$ shrinks faster than the sphaleron energy grows. When the Higgs mass 
reaches 6.665 Mw, the action at point Q crosses below the minimum action of the periodic 
instantons /, that is, An/aw- This situation is shown in figure |8|(c). 

Given this information available from perturbative calculations about the branches of 
periodic instantons and bounces, one can construct a likely scenario for the behavior of the 
various solutions in between the known limits. 

For a sufficiently light Higgs mass, where the situation shown in figure |8|(a) holds, it 
is obviously simplest to assume that the branch of periodic instanton solutions I and the 
branch of bounces C are one and the same. We illustrate this hypothesis in figure |9](a), 
where the dotted line indicates the single hypothetical branch of solutions. Both of the 
solutions / and C have two parity odd, time-reversal even negative modes, which is an 
important consistency check on this picture. This picture of a single branch of periodic 
solutions interpolating between the short-period instanton-antiinstanton solution and the 
sphaleron at its critical period is precisely the situation which occurs in the two-dimensional 
non- linear 0(3) sigma model with a symmetry breaking mass term ||. 

When the Higgs mass exceeds the value 3.091 M w , as discussed above, the bounce period 
increases as the turning point energy is decreased from the sphaleron energy. The simplest 
explanation for this, drawn from analogy with the toy model, is shown in figure ||(b). A 
new bifurcation point P appears, where the branch of bounces C and the branch of periodic 
instantons / merge, at a period longer than /3 - 

When the Higgs mass reaches the value 6.665 M w , as noted above, the action at point Q 
is equal to the minimum action of the periodic instantons /. This situation is illustrated in 
figure |9](c). This implies that, to exponential accuracy, the bounces set the rate of barrier 
penetration for the range of periods for which their action is below the periodic instanton 
action. 

As discussed earlier, bounces emerging from the sphaleron with f3 > (3q have only a 
single negative mode (which is parity odd and time-reversal even). Hence the branch QP 
for Mh > 3.091 Mw will have only a single negative mode. Point P is thus a typical example 
of a bifurcation where real solutions with stabilities differing by ±1 merge. 
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Period 
a) M H < 3.091 M w 



Period 

b) 3.091 M w < M H < 6.665 M w 




FIG. 9. The simplest scenario for Euclidean classical solutions of the Si7(2)-Higgs model consistent 
with the known perturbative results of figure |[ The number next to each branch of solutions gives the 
number of unstable directions, (a) For Higgs mass lighter than 3.091 Mw, the periodic instantons and 
bounces are joined by a single branch of solutions, indicated here by a dotted line, (b) For a Higgs mass 
above 3.091 Mw, the periodic instantons and bounces meet at a new bifurcation point P. (c) When the 
Higgs mass exceeds 6.665 Mw, the action at point Q is less than the action of the periodic instantons /. 

To exponential accuracy, the rate of barrier penetration at a given temperature in 
each of the plots in figure ^| is set by the lowest action configuration with period (3, except 
that there is an additional lower limit set by the zero size limit of periodic instanton configu- 
rations, which, for any period, have a limiting action of Aw/aw Thus, in cases (a) and (b), 
there is an abrupt crossover from a sphaleron-dominated rate at higher temperatures to 
a singular instanton-dominated rate at lower temperatures. In case (c), there is a smooth 
transition from the sphaleron to the branch QP, followed by an abrupt transition to singular 
instantons. 



NOTE ADDED 



The above scenario for periodic Euclidean solutions of S77(2)-Higgs theory has recently 
been confirmed by direct numerical studies of time-dependent (spherically symmetric) peri- 
odic solutions |0|[2(|. 
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